Evolution of entanglement entropy following a quantum quench: 
Analytic results for the XY chain in a transverse magnetic field 
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The non-equilibrium evolution of the block entanglement entropy is investigated in the XY chain 
in a transverse magnetic field after the Hamiltonian parameters are suddenly changed from and to 
arbitrary values. Using Toeplitz matrix representation and multidimensional phase methods, we 
provide analytic results for large blocks and for all times, showing explicitly the linear growth in 
time followed by saturation. The consequences of these analytic results are discussed and the effects 
of a finite block length is taken into account numerically. 

PACS numbers: 03.67.Mn, 02.30.Ik, 64.60.Ht 



The non-equilibrium evolution of extended quantum 
systems is one of the most challenging problems of con- 
temporary research in theoretical physics. The subject 
is in a renaissance era after the experimental reahzation 
[H of cold atomic systems that can evolve out of equilib- 
rium in the absence of any dissipation and with high de- 
gree of tunabihty of Hamiltonian parameters. A strongly 
limiting factor for a better understanding of these phe- 
nomena is the absence of effective numerical methods to 
simulate the dynamics of quantum systems. For meth- 
ods like time dependent density matrix renormalization 
group (tDMRG) this lack of efficiency has been traced 
back [3| to a too fast increasing of the entanglement en- 
tropy between parts of the whole system and the impos- 
sibility for a classical computer to store and manipulate 
such large amount of quantum information. 

This observation partially moved the interest from the 
study of local observables to the understanding of the 
evolution of the entanglement entropy and in particular 
to its growth with time [5|j6| . Based on early results from 
conformal field theory [6|,|7|]and on exact/numerical ones 
for simple solvable model @, Q it is widely accepted 
that the entanglement entropy grows linearly with time 
for a so called global quench (i.e. when the initial state 
differs globally from the ground state and the excess of 
energy is extensive), while at most logarithmically for a 
local one (i.e. when the the initial state has only a local 
difference with the ground state and so a little excess of 
energy). As a consequence a local quench is simulable by 
means of tDMRG, while a global one is not. 

However, despite this fundamental interest and a large 
effort of the community, still analytic results are lacking. 
In this letter we fill this gap providing the analytic ex- 
pression of the entanglement entropy at any time in the 
limit of a large block for the XY chain in a transverse 
magnetic field described by the Hamiltonian 
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where tr" are the Pauli matrices at the site j. Periodic 



boundary conditions are always imposed. In spite of its 
simplicity, the model shows a rich phase diagram being 
critical for /i = 1 and any 7 and for 7 = and \h\ < 1, 
with the two critical lines belonging to different univer- 
sality classes. The block entanglement entropy is defined 
as the Von Neumann entropy Si — —Trpglog pi, where 
pi = Tr„>^/9 is the reduced density matrix of the block 
formed by i contiguous spins. In the following we will 
consider the quench with parameters suddenly changed 
at time t = from /iq, 70 to /i, 7. 

Our main result is that, in the thermodynamic limit 
N 00 and subsequently in the limit of a large block 
£ 3> 1, the time dependence of Se{t) can be written as an 
integral over the momentum variable 
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where e' — de/dLp is the derivative of the dispersion re- 
lation = (/i — cos Lp)"^ + 7^ sin^ ip and represents the 
momentum dependent sound velocity (that because of 
locality has a maximum we indicate as vm = maxip |e'|), 
cos A;p = {hhQ — cos p{h + ft-o) + cos^ p + 770 sin^ p) / eeo 
contains all the quench information Q and H{x) = 
-[(1 + x)/2 log(l -f x)/2 + x)/21og(l - x)/2]. 

We first prove ([2]) and then discuss its interpretation 
and physical consequences. The readers not interested in 
the derivation can jump directly to the latter part. 

The method. The entanglement entropy can be writ- 
ten in terms of a block Toeplitz matrix [y, [l3| . One first 
introduce Majorana operators 02;-! = (nm<; '^m) "'f 
and 02/ = (nm<; '^m) "'f ^"^^ correlation matrix 
Vf through the relation {dm(in) = Smn + 
1 < m,n < £, that is a block Toeplitz matrix 
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cos6'(p = — (/i — cosiy9)/e, and sinAip = — sin (^[7/10 — 70/1 — 
cos (/3(7 — 7o)]/eeo- The entanglement entropy is given by 
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This trace can be numericaUy evaluated for finite H. as 
done in Ref. in the Ising case (7 = 70 = 1) for 
^,^0^1- A strong numerical evidence supports the 
fact that for 2vM't. < i, Se{t) increases linearly with time 
for large enough i, but this remained without any proof 
until now. It has also been argued that the limit t 00 
exists and Se saturates to a value proportional to £ (op- 
positely to the ground-state expectation where there is 
at most a logarithmic £ dependence 13, U, that has 
been calculated with a generalization of the Szego lemma 
. Eq. (I2|) not only provides the proof for a strictly lin- 
ear increasing of Si{t) for t < £/2vm, but gives also the 
complete time dependence. 

Proof of Eq. 0). Let us first sketch the strategy to 
prove Eq. ([2]) and give only after the technical details. 
The matrix iTg has the same eigenvalues of the £ x £ 
hermitian Hankel-|-Toeplitz matrices W± — H ± iT, with 
Hkj = ge+i-k-j and Tkj = fk-j- In fact, if w is an eigen- 
vector of W± with eigenvalue then the vector u with 
elements U2j-i = Wj and U2j — ±Wi+i-j is an eigenvec- 
tor of iTi with eigenvalue ^uj. We will show that Tr W^^ 
with n integer satisfies a relation similar to Eq. ([2|) with 
H{x) replaced by the appropriate power. Since Eq. ^ 
can be written as an expansion in Tr VF" with only even 
powers, this implies Eq. ([2]). Another ingredient is that 
in Eq. ^ for £ ^ 1 the term e~**»' is stationary being 
£ independent. At this point a multidimensional inte- 
gral for TrVF^" (see below) can be calculated with the 
stationary phase approximation, that is exact for £ ^ 1. 

In order to get Tr[T/F"]/^, n e N as ^ ^ cxd with a mul- 
tidimensional stationary phase approximation (see e.g. 
[3]), we consider the Toeplitz and the Hankel symbols 
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Each multiplication between two H - 
sums like 



' T matrices involves 
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where the sign is plus when an Hankel matrix is on the 
left and minus otherwise. We can fix the sign to be minus, 
multiplying any Hankel symbol to the right by the par- 
ity operator acting as Pf{x)P = f{—x). When we close 
the chain with the trace operator, we cannot change the 
sign of the first symbol and if there is an odd number 
of P, the term does not contribute to the leading order, 
because the phase in the integrand is not stationary any- 
more. The ^ dependence on the phase implies that all ipj 
variables are equal along the stationary curve. Note that 
all the significant terms have the same number of Hankel 
symbols with reversed signs and with not. In fact, con- 
sidering a tensor product of n symbols inside a symmetric 
integration, it follows (with ~ we always mean equal in 
the limit of large £) 

fc 
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with /i_ = P/iP, as can be straightforwardly proved 
by induction. Using the parity h*(p) = h{—p>), we 
have that the Hankel symbol phase e"*^"^, which is not 
proportional to £ can be dropped. The Toeplitz sym- 
bol is odd t{p) — ^t{—p) so the symmetrized prod- 
uct {t{pi) + h{pi)P) (g, {t{p2) + h{(p2)P) ~ t{(pi)t{ip2) + 

h(pi)h{—p2) leaves no parity terms, and if n is odd the 
whole integral asymptotically vanishes. Thus, from now 
on, we use the redefined symbols 
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which depend on the initial parameter only through A^. 
Repeated application of the multiplication rule leads to 
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where C is the hypercubic domain. The product A can 
be moved inside the cos function turning it in a sum, 
because of the symmetry of the ^ domain of integration 
with respect to 0. There is a trivial integration along 
a direction in the ^ domain since the integrand depends 
only on the difference between the ^ variables, thus 
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= max[0, min(l, 1 - V Ci) + min(l, 1 + V CO] • 

The permutation symmetry of ip variables allows to order 
the cos products in _B, so that we can introduce a set of 
spin variable aj G { — 1,1} and bring B to the form 
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Interchanging the limit £ — ^ oo with the integration with 
respect to the variable ip = ^p2m the remaining (2n — 2)- 
dimensional integral is easily solved by stationary phase 
methods. For each configuration of {a} the Hessian de- 
terminant is 2^^"'" and the Hessian signature vanishes 
(this is the reason why in Eq. ([2]) oscillations, usually 
present in stationary phase calculations, are not present). 
On stationary points we finally have 11 = 2 on time- 
independent terms and Jl = 2(1 — 2\e'\t/t) otherwise. Of 
course the limit exists and the direct computation gives 
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that is a "Taylor" expansion of Eq. Q . From Eq. ([1 
we also have the time dependence of all Renyi entropies 
Sr = logTr[p"]/(l — a): it is enough to replace i/(cos A) 
with logdcos^ f I" -t- Isin^ A|")/(i _Q,) in Eq. ^. On 

passing, we mention that our exact result obviously sat- 
isfies the bound for Si(t) given in Ref. [1]. 

Description of the result. In Ref. an interpretation 
of the time dependence of Si has been provided in terms 
of causality (later generalized to the correlation functions 
in jJL3i])- The idea is simple: the initial state has a very 
high energy relative to the ground state of the Hamilto- 
nian which governs the time evolution, and therefore acts 
as a source of quasiparticle excitations. Particles emitted 
from different points (further apart than the correlation 
length in the initial state) are incoherent, but pairs of 
particles moving to the left or right from a given point 
are entangled. Thus Si{t) should just be proportional to 
the number of coherent particles that emitted from any 
point reach one a point in [0, t] and the other the remain- 
der of the system. Since there is a maximum speed for 
these excitations vm, this implies the linear growth for 
2vMt < i and saturation for very large times. 

However, only in the conformal case when e' does not 
depend on the momentum because of the linear disper- 
sion relation, this scenario makes quantitative predictions 




FIG. 1; R defined in Eq. ((13} for quenches from {hit = 
10, 7o = 2) . The inset shows the singular behavior in the 
neighborhood of /i = 1. 



on the time evolution, else the rate of production of parti- 
cles f{p',p") is an unknown function of the Hamiltonian 
parameters both before and after the quench. The com- 
parison of Eq. (21) with the general one (Eq. (4.2) in Q) 
allows to identify /(p',p") with i5(p' -p")iJ(cos Ap/). We 
can also easily read from our result the value of the ratio 
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that results to be the average of the absolute value of 
speed of the sound on the H(cos A^) distribution. R as 
function of the quench parameters is shown in Fig. [TJ It 
is not analytic at the quantum critical point /i = 1, as a 
trivial consequence of the non-analyticity of its building 
blocks (i.e. e, A^). However it is clear from the inset 
that such non-analyticity is so weak that is unrealistic to 
say that the out-of-equilibrium behavior of entanglement 
entropy is sensitive to the phase transition. 

From Eq. ^ we also have the large time corrections 
to the asymptotic result. Since H{±.1) = with a log 
singularity, when the zero- velocity mode giving the large 
t behavior is at = ±7r (as e.g. for h > 1), one has that 
the first correction is oc log t/t^, whereas when there 
are zero- velocities not at the border of the Brillouin zone, 
where H{x) is finite, the leading correction is oc £^/t. 

For t = oo only the second term in Eq. ^ contributes 
to the entropy that thus is extensive. As already no- 
ticed in Ref. Q this is the same result at finite large 
temperature PeS- An interesting question is whether this 
effective temperature is observable independent as found 
in the conformal case [l3| or instead depends on the oper- 
ator and so would not be a well-defined concept. Further 
checks of this point are mandatory to avoid speculations. 

As a consequence of H{x) < H{0) ~ In 2, we have 
S£{t) < £\n2 for any time, a bound that is just the 
maximum entanglement allowed by the dimension of the 
Hilbert space. In Ref. [gJ it was noticed that the various 
curves for the quench from /iq = oo to any h apparently 




FIG. 2: Time evolution of the entanglement entropy Si{t)/l 
for several quenches and I. The straight line is the leading 
asymptotic result for large I. The inset in the bottom-left 
graph shows the derivative with respect to time of Se{t) for 
I ^ oo and the numerical derivative for I = 90. 

collapse on a single curve when rescaling Si{t) to ^'^(oo). 
From Eq. this is exactly true for \h\ < 1, but only 
approximately otherwise. Finally, the t — oo result is 
symmetric under the exchange {h, 7) (/iq, 70), because 
the asymptotic result only depends on cos A^p that does 
not distinguish between initial and final values. 

Finite £. The matrix representation Q allows for 
the numerical calculation of Se for finite and relatively 
large £. Some results are reported in Fig. ^ where we 
choose those quench parameters that make the finite £ 
effects more relevant. Increasing i the results always ap- 
proach Eq. showing unambiguously its correctness, 
but there are peculiar and interesting finite £ effects. The 
most evident effect is the oscillation of Si{t). These os- 
cillations have been generically seen in numerical studies 
also on more complicated models P] , but they are absent 
in the Ising model for h,ho > 1 [G]: oscillations can only 
be present when there is a second local maximum of |e'|. 
The data (also for cases not shown in the figure) provide 
a strong evidence that the first non-oscillating correction 
at order 0{£°) is positive and time independent. 

In the bottom- left plot in Fig. [5] the most unexpected 
effect is shown. For the quench {Iiq — 1,70 = 0.4) — > 
{h — 0.5,7 — 0-2), it seems that the linear regime of 
Si{t) continues after t* = £/2vm- However, looking at 
the derivative (inset) one realizes that it is not exactly 
constant, since it slightly bends at t* . This happens 
because for this peculiar quench the maximum velocity 
mode carries very little information, and so a stronger 
non-analyticity is present at a local maximum of the ve- 
locity smaller than vm- This effect is pronounced every 
time that hoj ~ /170, with \h\, I7I, |/io|, I70I < 1, because 
of the functional form of A^. This anomalous behavior 
is important because it is nowadays common to extract 



the speed of propagation of information from t* . Every 
time this effect is present, this procedure gives the wrong 
answer. For example, we plotted in the inset of Fig. [2]the 
numerical derivative of Se{t) for £ — 90 {a, value hardly 
reached in non-equilibrium simulation). It is evident that 
at t* there is no trace of the non-analyticity. Relying on 
these results one would have obtained a value of vm that 
is almost half of the real one. 

Conclusions. The non equilibrium time-evolution of 
Se(t) for the XY chain seems to encode most of the fea- 
tures that have been observed numerically in other con- 
texts It is then natural to wonder whether slight 
modifications of Eq. ^ can be true in more compli- 
cated situations and not only for models mapped to free 
fermionic theories as the present one. 
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